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We show how particle-vortex duality implies the existence of a large non-abelian discrete symmetry 
group which relates the electromagnetic response for dual two-dimensional systems in a magnetic 
field. For conductors with charge carriers satisfying Fermi statistics (or those related to fermions by 
the action of the group), the resulting group is known to imply many, if not all, of the remarkable 
features of Quantum Hall systems. For conductors with boson charge carriers (modulo group trans- 
formations) a different group is predicted, implying equally striking implications for the conductivi- 
ties of these systems, including a super-universality of the critical exponents for conductor/insulator 
and superconductor/insulator transitions in two dimensions and a hierarchical structure, analogous 
to that of the quantum Hall effect but different in its details. Our derivation shows how this sym- 
metry emerges at low energies, depending only weakly on the details of dynamics of the underlying 
systems. 
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I. INTRODUCTION 

Two-dimensional electron systems have remarkable 
properties, including the quantum Hall effect and metal- 
insulator transitions, many features of which still resist 
theoretical explanation. The difficult part of describing 
these systems is that they involve strong correlations, and 
no small parameters present themselves to help with the 
analysis. This deprives theorists of most of the tools in 
their conceptual toolboxes. 

Two kinds of theoretical tools which have proven use- 
ful for analysing these kinds of strongly- coupled problems 
are the exploitation of symmetries and of the simplifi- 
cations which are associated with the low-energy limit. 
Duality symmetries, in particular, are likely to be useful 
since these typically relate strongly-coupled degrees of 
freedom to weakly-coupled ones, and in two dimensions 
particles and vortices make natural candidates for dual 
partners. Indeed, particle-vortex duality has been used 
several times in the literature to describe some aspects 
of both the Quantum Hall effect and conductor-insulator 
transitions in superconducting films. 

Duality symmetries arc unusual in that they are not 
symmetries in the usual sense that they need commute 
with the system's Hamiltonian. Instead, they relate two 
different kinds of systems to one another. This relation- 
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ship is useful when one of the two systems so related can 
be analysed, permitting conclusions to be drawn for its 
dual by acting with the duality transformation. Alterna- 
tively, duality can be useful if it maps a family of systems 
into themselves, since duality invariance then constrains 
how one flows amongst members of the family as exter- 
nal parameters - like temperature or magnetic field - are 
varied. 

Our aim is to show that both of these lines of argu- 
ment have very general applications to two-dimensional 
systems. In particular, our main result is to show that 
the twin operations of particle-vortex duality and 
the addition of 27r units of statistics to charge carriers 
(which does not change their statistics at all) take a very 
simple form when expressed in terms of the material's 
electromagnetic (EM) response functions.^] These rela- 
tions hold for any system for which the low-energy EM 
response is dominated by the motion of quasi-particles 
or vortices, and for which the dynamics of these quasi- 
particles and vortices are similar (in a way we make more 
explicit in what follows) . Since these two transformations 
do not commute, they generate an infinite discrete group 
of duality relations amongst the EM response functions, 
and it is this large group which underlies the predictions 



* Systems for which dualities commute with the Hamiltonian 



make up the special case of self-dual systems. 

^See references |^-^| for other approaches to duality in Quan- 
tum Hall systems. 
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which we shall find. 

In general the duality transformations we find do not 
preserve the momentum dependence of the EM response 
function and so, for instance, can relate materials whose 
response is very different (such as by relating supercon- 
ductors to insulators). It turns out that conductors are 
mapped into themselves, however. When specialized to 
conductors, with the electromagnetic response character- 
ized by the Ohmic and Hall conductivities, a xx and a xy , 
the action of duality takes the form of subgroups of the 
modular group, PSL{2, Z), acting on the complex con- 
ductivity, a = <r xy + icr xx f\ as follows: 



with the integers a through d satisfying ad — be = 1 (the 
occurrence of this symmetry in a statistical mechanical 
model was first noticed in ||, in an investigation aimed 
at understanding QCD). The duality transformations as 
defined in eq. (|l|) are not symmetries of the Hamiltonian 
since, for instance, dual pairs differ in their electromag- 
netic response. Rather, these transformations are sym- 
metries of the flow obtained as external parameters, such 
as magnetic fields and temperatures, are varied. 

Which particular subgroup is important depends on 
the statistics of the charge carrying quasi-particles of the 
problem. If they are fermions, or related to them by 
the symmetries we shall describe, then the subgroup (de- 
noted To (2) in the mathematics literature jl(i|0 is defined 
by the condition that the parameter c must be even. For 
bose charge carriers, or their symmetry partners, such 
as for superconducting films or Josephson Junction ar- 
rays, the symmetry is instead 1^(2), defined by the con- 
dition that a, d are both odd and b, c are both even, or 
vice versa. For both cases powerful predictions follow 
from the fact that our derivation shows that the symme- 
try constrains how the conductivities change as external 
variables are varied. 

For Quantum Hall systems, the group To (2) has been 
long conjectured to be important and has been 

derived for these systems within a mean-field approxima- 
tion For these systems our derivation accomplishes 
two new things. First, although our arguments are mod- 
elled on those of ref. [Q, ours have a broader domain 
of validity since they explicitly assume only that quasi- 
particles or vortices dominate the low-energy EM re- 
sponse, and that the systems are clean enough to exclude 
any interactions which might distinguish quasi-particles 
from vortices, and so thereby ruin the duality symmetry 
which relates them. Because the fields we treat only arise 
in the effective theory, and are not meant to describe the 
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complete electron dynamics, the mean-field approxima- 
tion is kept under better control. 

Second, our derivation helps clarify the assumptions 
which underlie analyses of the consequences of Fq(2) 
invariance for the renormalization-group (RG) flow in 
the o~ xx — <j xy plane, since we show that this only re- 
lies on the underlying particle-vortex duality and on the 
long- wavelength limit. This is important because it has 
been shown [p^|-p^[ that most of the unique features 
of Quantum Hall electromagnetic response follow from 
the consistency of To (2) invariance with RG flow in the 
&xx — o~xy plane, independent of the detailed form of 
the flow's /^-function. (The constraints on (3 which fol- 
low from this symmetry have also been considerably ex- 
plored Jl4|,[l6| |ij| .) Previously the key assumption of two- 
dimensional flow, governed by To (2) invariance, was just 
that: an assumption, although a plausible one motivated 
by analogy with the two-dimensional scaling theory of 
disorder [po|-p3|. In particular, since the scaling theory 
strictly only applies near the system's critical points the 
scaling motivation could not explain why many of the 
predictions following from ro(2)-invariant flow work ex- 
tremely well, even away from the flow's critical points. 
Since the derivation presented here is not similarly re- 
stricted to scaling regions, it explains why these other- 
wise surprisingly successful predictions work. 

The identification of r#(2) as the duality symmetry for 
two-dimensional conductors with bosonic charge carriers 
is new, although some precursors of this idea exist. Im- 
plications of the particle-vortex generator of the group 
for critical behaviour in a superconductor-insulator tran- 
sition have been examined pi, and the group r#(2) was 
discussed as potentially playing a role in two-dimensional 
systems |ll]|n| , including possibly for the Quantum Hall 
effect ||. The action of Tg(2) on longitudinal conduc- 
tivities was written down in p4| , although this was not 
extended to the whole upper-half complex plane. 

Also new to this paper are the detailed predictions 
which follow for bosonic systems from the proposed sym- 
metry group, and which arc direct analogs of the symme- 
try consequences which are already known for Quantum 
Hall systems. 

To which real-life systems does our duality symmetry 
apply? Just as it is difficult to compute ab initio whether 
a material will be a solid or not, it is difficult to an- 
swer from first principles which systems must be particle- 
vortex symmetric. There are several things which can be 
said, however. First, if the symmetry holds, then all of its 
consequences follow together. For example, for Quantum 
Hall systems the semi-circle law, p xx — > l/p xx duality, 
super-universality of the critical exponents, universality 
of the transition conductivities, odd-denominator frac- 
tional quantization of o~ xy when o~ xx — 0, the selection 
rules for which plateaux may be related by transitions, 
etc. must all come together as a package. 

Second, since particle-vortex duality relies on the 
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equivalence of the kinematics and couplings of the charge- 
carrying quasi-particles and vortices, it should be a good 
approximation when the only quantities of interest in the 
Hamiltonian are those describing the kinematics of these 
particles, and their couplings to the fields which describe 
the long-range vortex interactions and the electromag- 
netic fields which are applied to test the EM response. 
Duality could be ruined by other microscopic interactions 
which treat quasi-particles and vortices differently, such 
as from couplings with disorder or with other electronic 
degrees of freedom. Of course, disorder can also play 
other spoiling roles, such as by destroying the phase co- 
herence on which the quantum regime which we assume 
depends. 

A sufficient condition for particle-vortex duality, and 
its associated non-abelian extensions, might therefore be 
that the system be sufficiently clean to justify the neglect 
of other particle and vortex interactions when calculating 
the electromagnetic response. Although this condition is 
not strictly necessary, since our derivation also applies 
in the presence of any particle/ vortex-democratic inter- 
actions, weak coupling is also implicit in our neglect of 
anomalous dimensions when deciding the relevance or ir- 
relevance of low-energy interactions. Notice that it is the 
weakness of the quasi-particle couplings which are impor- 
tant in this decision, and the assumption that these are 
weak does not imply that the underlying electrons must 
also be weakly coupled in the microscopic Hamiltonian. 

Now to the main arguments. We organize our presen- 
tation as follows: First we describe the action describing 
the low-energy dynamics of a system of quasi-particles 
and vortices, and cast it into a form which emphasizes 
the similarities between these two kinds of charge carri- 
ers. Next, we derive the action of the two basic symme- 
tries — 2ir statistics addition and particle- vortex duality 
— for the electromagnetic response functions. Then, we 
specialize the result to the particularly interesting case of 
a conductor, to derive the action of the symmetry on the 
conductivities, o~ xx and o~ xy . We briefly review the Quan- 
tum Hall case, where the charge carriers are fermions, 
and then repeat the analysis for charge carriers satisfy- 
ing bose statistics, listing many experimental predictions 
which follow from the symmetries. 

II. PARTICLES AND VORTICES 

This section has two goals. First, we derive an ex- 
pression for the effective action governing the low-energy 
interactions of charged quasi-particles and vortices for 
which the duality between these two kinds of objects is 
made explicit. Second, we compute the EM response for 
a system of such charges and vortices, for use when de- 
riving the implications of particle-vortex duality in the 
next section. 



A. The Effective Quasi- particle- Vortex Action 

Our starting point is the following lagrangian, which 
describes the low-energy/long- wavelength interactions of 
a collection of N p quasi-particles and N v vortices with a 
weak electromagnetic probe, A^: 

£eff = - ^ ^ Xv a^d x a, + £ p (£, a + A) 

+ •••, (2) 

which we write in the continuum approximation because 
our interest is directed towards the low-frequency, long- 
wavelength EM response. 

Several features of this lagrangian bear explanation, 
and since our final results ultimately depend on its va- 
lidity, we pause here to explain its form in some detail. 

• Why Both Particles and Vortices?: Usually charged 
particles and vortices, in the way we define them, do 
not coexist in the low-energy theory, since the vor- 
tices presuppose the breaking of electromagnetic symme- 
tries which precludes the existence of isolated electrically 
charged particles, and our later applications only require 
the consideration of systems containing one or the other. 
We nonetheless use the mathematical device of keeping 
both in eq. (||) since it permits us to derive our results for 
particles only and vortices only by taking the appropriate 
limits of a single formula. 

• Field Content: A^ is the small external electromagnetic 
potential which is applied to probe the system's EM re- 
sponse, £fcW is the position of the ft'th quasi-particlc of 
the system as a function of time, and a M is the usual 
statistical gauge potential which ensures that the inter- 
change of two quasi-particles produces the phase e t0 [^5) . 
(9 = 2nn, for integer n, corresponds to bosonic quasi- 
particles, while 9 = (2n + l)ir describes fermionic quasi- 
particles. |] 

• Relativistic Form: We use relativistic notation in eq. (|^) 
in order to most cleanly illustrate the logic of the argu- 
ment. This proves to be convenient because the relativis- 
tic case shows all of the main features of duality, and is 
considerably simpler to describe. We have checked that 
other features of a non-relativistic treatment - such as 
the potential appearance of kinetic terms linear in time 
derivatives - do not substantially affect our arguments. 
It also happens that the results for the non-relativistic 
systems of practical interest can be read off directly from 



^Appendix A briefly reviews our conventions concerning this 
statistics field. 
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the relativistic answers using the trick outlined in Ap- 
pendix B. 

• Quasi-particle Lagrangian: £ p (£,a + A) is the (first- 
quantized) lagrangian which describes the quasi-particle 
motion and their coupling to the electromagnetic field, 
A^. A first-quantized representation is chosen because 
this makes the duality between particles and vortices 
most transparent later on. The detailed form of this la- 
grangian is not important in what follows, apart from the 
form of the coupling to the gauge field, (a + A)^, whose 
form is important, but which follows on grounds of gauge 
invariance. To be concrete, in the absence of other inter- 
actions the particle lagrangian might be given explicitly 
by 



[ 2 £fc£ fc A< 



qk&ia + A)^ 5[x-£k(t)], 



(3) 



where qk here denotes the quasi-particle charge, normal- 
ized so that qk = — 1 for electrons. 

• Dependence on External Parameters: All of the depen- 
dence on the external variables, like magnetic field, B, 
enter through the parameters of the effective lagrangian. 
For instance, in the example where the quasi-particle dy- 
namics is described by eq. (||) they would enter through 
the quasi-particle mass, m, and the parameter n, al- 
though the dependence would be more involved for more 
complicated quasi-particle dynamics. The same is almost 
true for the dependence on temperature, T. That is, 
the temperature dependence contributed when integrat- 
ing out high-energy modes is embedded in the system pa- 
rameters, but there is also additional temperature depen- 
dence associated with the integration over the low-energy 
degrees of freedom themselves, such as Both forms of 
temperature dependence are included in our later discus- 
sions of the temperature-dependence of particle-vortex 
duality. 

• The Goldstone Variable: If the difference between the 
number of vortices and anti- vortices, N v , is nonzero then 
there is a complex order parameter somewhere in the 
system which typically vanishes at the positions of the 
vortices (and anti- vortices) , and takes a nonzero value 
asymptotically far away. The number N v is then re- 
lated to the winding of the phase of this order parameter 
around a circle which encloses all of the vortex positions. 
The field cj> in eq. (|J) represents the phase of this order 
parameter. We assume, in writing eq. (^|), that the order 
parameter carries nonzero electric charge, q^ ^ 0, and so 
it spontaneously breaks the electromagnetic U(l) gauge 
group. The scale of the parameter n is of order of the 
scale of the symmetry-breaking expectation value. 

Since the quanta of <j> are the Goldstone bosons for 
the assumed symmetry breaking, <f> is guaranteed to be 
in the low-energy theory. Indeed <p would be responsi- 
ble for the long-range interactions experienced between 



vortices, if q^ were zero. Furthermore, its couplings to 
(a + A)^ are dictated by gauge invariance. Finally, be- 
cause it is a Goldstone variable 4> is guaranteed by Gold- 
stone's theorem to decouple at low energies (modulo the 
usual Coleman-Mermin- Wagner caveats), thereby justi- 
fying its semi-classical treatment using a derivative ex- 
pansion |^(| . The ellipses in eq. (^) represent all of the 
other effective interactions obtained when all higher de- 
grees of freedom are integrated out. Since these all in- 
volve inverse powers of the higher-energy scales, such 
as the order parameter scale, they are irrelevant for the 
present purposes compared to those explicitly displayed. 

In the derivation of KLZ [fijl this order parameter was 
the bosonic field which described the electrons, but for 
the present purposes it could equally well be a bilinear 
of fermion fields, or something still more complicated. 
All we need assume is that the order parameter spon- 
taneously breaks electromagnetic gauge invariance, and 
that its boundary conditions at spatial infinity incorpo- 
rate the winding corresponding to vortex number N v . 



B. A Dual Description of the Vortices 

In order to better display the particle-vortex duality, 
it is convenient to use the dual description of the vor- 
tex degrees of freedom |^-^]. We here implement this 
duality transformation by recognizing that it is a special 
case of a general dualization algorithm [B(J (which also 
has applications to bosonization in one and two dimen- 
sions and is derived in detail for the present system 
in Appendix C. (We repeat the derivation in a second- 
quantized format in Appendix D, using the same model 
- the nonrelativistic abelian Higgs model supplemented 
with Chern Simons coupings - used by KLZ.) 

The result, when applied to the quasi-particle/ vortex 
action, eq. (0), is: 



a + A) + £ v (y,b) + . . . , 



(4) 



where: 



• The field 6 M : is the new (gauge potential) field which 
is the dual representation of the Goldstone field, (f>. It 
carries all of the information about the long-range inter- 
actions amongst the vortices. 

• The variables y% (t) : label the positions of the centres of 
the vortices and anti- vortices defined, for instance, as the 
positions of the zeros of the underlying order parameter 
whose phase was </>. As is seen in Appendix C, these 
positions naturally arise as variables during the duality 
transformation once one takes into account the boundary 
conditions satisfied by <fi in the presence of vortices. 

• The Vortex Lagrangian C v (y,b): describes the dynam- 
ics of the vortices and their couplings to the field o^. 
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Although this action can be complicated, reflecting the 
potentially complicated dynamics of vortices in the ma- 
terial of interest, it must have at least the following two 
terms: 



2 VaUafl 



2nN a 



6[x-y a (t)} ) (5) 



where N a denotes the vorticity (or winding number) of 
each vortex. The coupling term between y£ and 6 M falls 
directly out of the duality transformation, and so is quite 
generally known. It should be noticed that although the 
kinetic term for y% does not itself follow directly by dual- 
izing eq. (^|) , its form is fixed quite generally from symme- 
try arguments ]3^ |. This is because the y£"s may them- 
selves be thought of as Goldstone bosons for the breaking 
of translation invariance, which is here broken by the po- 
sitions of the vortices themselves. In general there may 
also be other interactions to supplement eq. which 
describe the interactions of the vortices with other de- 
grees of freedom, such as disorder. 

A b u kinetic term: of the form -Trbr ft. v {f b Y v ■* witn 



fuu = d^by — dvbft, is also produced when perform- 
ing the dualization, but is not written in eq. (Q). It 
has been dropped since it is inversely proportional to 
k, which is one of the high-energy scales whose in- 
verse we are ignoring in the low-energy, long-wavelength 
limit. It would compete with a term proportional to 
(fnv + F^)(P V + Fn (where and are the field 
strengths for and A^), which was among the ellipses 
appearing in eq. (^) . (Alternatively, both such terms may 
be absorbed into the general response function which is 
obtained when and y£ are integrated out, as we are 
shortly to describe.) 

• Equivalence: It is to be emphasized that eq. (B) is just 
a change of variables of eq. (0) , and so describes precisely 
the same physics. In particular, both expressions repro- 
duce precisely the same electromagnetic response once all 
degrees of freedom except for A^ are integrated out:|**| 



[<M*)] JJ[de£(i)] exp i J d 3 xC cS (t a, A) 

k L 

[da,(x)][db^x)] l[[dv£(t)] 

a 

i / d 3 xC c g(y,a,b, A) 



x exp 



(6) 



** Although we present our results for real time and at zero 
temperature, our duality relations can be derived equally well 
in imaginary time at nonzero temperature. 



III. PARTICLE- VORTEX DUALITY 

The beauty of eqs. (||), (Q) and (||) is that they dis- 
play the quasi-particle and vortex degrees of freedom in 
a way which emphasizes the similarity of the particles 
and vortices. For instance, if other interactions are neg- 
ligible, so the particle and vortex dynamics is given by 
eqs. (||) and (Q), and in the absence of the Chern-Simons 
term involving (9, eq. has the symmetry with respect 



to the interchanges «-> y%, b p 



Au, so long 



as the quasi-particle masses and charges are also inter- 
changed: qk «-> 2TrN a /q < f, and m «-> M. The same is also 
true if C p and C v are more complicated, provided that 
the additional complications treat particles and vortices 
democratically, by contributing terms of the same form 
to both C p and C v . 

We may now be more precise as to what is meant by 
dual systems. Given a two-dimensional system whose 
EM response is governed by n quasi-particles (having 
mass m etc) and N vortices (having mass M etc) , we de- 
fine the dual to be the system having N quasi-particles 
(with mass M etc) and n vortices (of mass m etc). The 
cases of real interest for the applications which follow are 
the cases involving particles only or vortices only: n = 
and N = 0. 



A. Expressions for the EM Response 

We now wish to determine how the electromagnetic re- 
sponse of dual systems are related to one another. To do 
so, imagine integrating out the quasi-particles and vor- 
tices by performing the path integral over their positions, 
and y^. If we are only interested in linear response, 
we need not actually evaluate these integrals, but may 
parameterize them in terms of response functions, 
and V*" , as follows. That is, if we define: 



J~f p (a+A) 



J IlKWl ex P 1 J d 3 x£ p (t,a + A) 

k 

= [ Y[[dv£(t)] exp i f d 3 xC v (y,b) , 
then for linear response it suffices to take: 
7p (a) =-- f d 3 xd 3 x' a^{x)P' tu {x - x')a v {x'), 



(7) 
(8) 

(9) 



lv{b) = --j d 3 xd 3 x' bp{x)V^{x-x')b u (x'), (10) 

where P liV {x — x') and V^ v {x — x') define the particle and 
vortex response functions. Notice that for dual systems 
we expect P" y = V» v and V y = P^ y ', where the tilde 
denotes the result evaluated in the dual system. 
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If eqs. (0) - (10) are used in (||), then the remaining 
integrals over a M and b^ are Gaussian and so may be eval- 
uated explicitly to obtain the electromagnetic response 
function, defined by: 



Y(A) 



d 3 xd 3 x' A^{x)W{x - x')A v {x'). (11) 



One might be queasy about the consistency of first ex- 
panding to quadratic order in & M and then integrating 
over all values, and this queasiness would be justified if 
the expansion to quadratic order was done because b^ is 
small. Such a calculation really presumes an effective- 
field theory approach, where all terms involving higher 
powers of fields are irrelevant (in the RG sense), and so 
may be neglected within the very low energy effective 
theory. This is true in the present instance so long as 
anomalous dimensions are small, so that relevance may 
be judged purely using naive dimensional analysis. It is 
here that we implicitly assume the quasi-particles and 
vortices to be weakly coupled to other degrees of free- 
dom, although this does not also imply weak coupling 
for the underlying electrons. In any case, we later pro- 
vide an alternative derivation of a subset of our results 
which does not rely on the quadratic approximation of 
eqs. @, ©. 

We believe that our treatment of the quasi-particles 
purely within an effective-theory framework represents 
an important conceptual difference between the argu- 
ments formulated here, and those presented in the spirit 
of KLZ, in Appendix D. They differ because in the KLZ 
framework the bosonic field whose phase describes the 
vortices is taken to be a direct description of the under- 
lying electrons, rather than a low-energy effective field. 
As a result a weak-coupling mean-field analysis for the 
KLZ field is more directly tied to the strength of the 
couplings of the underlying electrons. 

Before quoting the result obtained in this way for the 
EM response, it is worth first defining some notation. 
The polarization tensor is usually taken to have the stan- 
dard rotationally-invariant and gauge-invariant but non- 
relativistic form: 



T(A) 



1 



d 3 xd 3 x' 



TT^x - x')E z (x)Ei(x') 



-tt 2 (x- x')B(x)B(x) (12) 
-ir 3 (x-x')e^ x A^(x)d u A x (x'), 

defining the electromagnetic form factors, 7Ti,7T2 and 7r 3 . 
In what follows our main interest is in m and 7r 3 , which 
control the conductivities o~ xx and o~ xy . Because of this, 
and because of the greater simplicity of the resulting for- 
mulae, we specialize instead to the relativistic version of 
eq. ([12[), which we write in momentum space as: 



IP" = U^p 2 ) A^ + n 3 (p 2 ) J"", 



(13) 



where A^ ; 
and J» v -- 



and J^ v are defined by A M „ — 5£ — p^p v /p 2 



Pxh 



Because the form factors of 



eqs. ( |l2] ) and ([l3|) are related (in momentum space) by 

IT = {p/h^m = (p/h) 2 Tr 2 and n 3 = ^Jp 2 /h 2 tt 3 , the 
relativistic form is sufficient to follow how the quanti- 
ties 7Ti and 7T 3 transform under duality transformations. 
(More generally, results for the generic case wi ^ 7T2, are 
easily obtained using the trick described in Appendix B.) 

The great utility of the relativistic expression, eq. (|l3|), 
follows because the tensors h.^ v and J^ v satisfy the iden- 
tities: A M a A Q „ = A^„, h^ a J av = J^ a A» a = J» u and 
J m Jav = — A M „. Since the tensors A and J are related 
to one another in the same way as are the bases, 1 and 
i, of complex numbers, tensor manipulations with IP" 
can be greatly simplified by re-expressing it as a complex 
variable: 



n = it + m 3 



(14) 



With this notation, and defining similar expressions 
for the complex quantities P and V in terms of the cor- 
responding form factors Pi, P 3 , V\ and Vs, we find the 
result of integrating the fields a M and 6 M out of eq. (Q) to 
be 



n 



i \j p 



p 



VP 



p 



VP 



(15) 



as the general relation between the EM response and the 
particle and vortex response functions. To avoid explicit 
factors of fi we have defined p — p/h in this expression. 
(The details of these integrations are given in Appendix 
E.) This expression is the main result on which our later 
conclusions are based. 

We now record the special cases of this expression 
which will be used in the following sections. 

• Addition of2ir Flux: Since 9 represents the statistics of 
the quasi-particles, the choice 9 = 2fc7r, with k an integer, 
can have no physical effects: H(9 + 2kn) « n(#). Inspec- 
tion of eq. ( [l5|) shows the implications of this statement 
for the EM response: 



1 



1 



1 



2ki 



11(6) TL{9 + 2kir) 11(9) 



(16) 



which reproduces a well-known result fllz 



• Quasi-particles Only: Should there be no vortices par- 
ticipating in the EM response at all (which, because of 
the duality transformation from cj) to b tll is equivalent to 
V — > oo), then eq. (fL5^ simplifies to: 



n 



G) 



(17) 



which reduces to II = P when 9^0, corresponding 
to bose statistics for the charge-carrying quasi-particles. 
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The corresponding result for fermionic charge carriers is 
similarly found by choosing 9 = ±7r. 

• Vortices Only: In the event that only vortices are in- 
volved in charge transport, expression ( [l5|) reduces to: 

n= ,„, . (is) 



V-i- 



'V 2 I 



This also reduces to the usual result, II 
bosonic charge carriers, for which 9 — > 0. 



p 2 /V, for 



IV. SOME CONSEQUENCES OF 
PARTICLE- VORTEX DUALITY 

In this section we derive the implications of particle- 
vortex duality for the electromagnetic response of two- 
dimensional systems. 

We are interested in the situation where quasi-particles 
only or vortices only are responsible for charge transport, 
in which case a very simple expression may be derived for 
the action of particle- vortex duality on the EM response 
function, II. This may be derived by using eq. ( fT^ ) to 
relate II to P for the original system, and using eq. JT§| ) 
to relate II to V = P for the dual system. Eliminating 
P from these two expressions gives the desired relation 



rw (f)-n 



p 



p 



^(f + #) n 



(19) 



directly expressing the dual response function, II, in 
terms of that of the original system. 

The physical interpretation of expression ([19]) is ob- 
tained by using the connection between the small-p 2 limit 
of the form factors, II i and II3, and measurable electro- 
magnetic quantities. There are two main cases to con- 
sider: 



A. Superconductors and Insulators 



1. Superconductors 

If A\ is nonzero, then the material is a superconductor 
and A\ — m 2 — 1/A is the corresponding photon squared 
mass, or inverse penetration length. In this case there is 
a gap because of the Anderson-Higgs mechanism. 

Inserting assumption ( p0| ) for II (with A\ 7^ 0) into 
the duality expression, eq. (|l9|), implies H also has an 
expansion of the form (|20|), with: 



(21) 



Ai = 0, 
A 3 = T 

Bi = — 



1 



We see that the system dual to a superconductor is an 
insulator, whose Hall conductivity and dielectric function 
are related to the superconducting penetration length 
and the statistics of its quasi-particle charge carriers. 

• Bosonic Charge Carriers: A particularly interesting 
special instance of expression ( |2l| ) is the case of bosonic 
charge carriers (such as Cooper pairs), for which 9 = 0. 
In this case we find as dual a dielectric with vanish- 
ing Hall conductivity and dielectric function given by 
Bi = 1/A±. This is as expected on physical grounds 
|l| since the condensation of vortices in the dual system 
should produce an insulator. Notice that eq. ( |2l| ) makes 
specific predictions as to the dependence of the dual di- 
electric function, Bi, as a function of temperature since 
it is related to the temperature dependence of the pene- 
tration length in the superconducting system. 

• Fermionic Charge Carriers: If we instead choose 9 = 
7r, as is appropriate for fermion charge carriers, we find 
the dual system has no Ohmic resistance, but has Hall 
conductivity a xy = h. 



If there there is a gap at the relevant part of the spec- 
trum, then one expects the EM response function to be 
analytic in momentum space. This implies the functions 
Hi and n3 have the following small-p 2 form; 



IIi = A t +B 1 p 2 



n. 



P 2 (A 3 



(20) 



The 



does not conflict with Ti^ v being analytic as 



p M — > because it cancels a similar factor in the def- 
inition of the tensor J^ v . The constant A\ is nonzero 
only for superconductors, with A\ inversely related to the 
medium's electromagnetic penetration depth. If Ay = 
then the material is an insulator and B\ is related to 
the dielectric response of the medium. The constant A 3 
corresponds to a Hall conductivity for the system, which 
may vanish but need not. 

There are two important sub-cases to consider: 



2. Insulators 

If A\ =0 then the material is an insulator and the 
gap exists because of the absence of low-energy charged 
quasi-particles which carry current. In this case B\ is 
related to the dielectric response of the medium and the 
duality transformation produces the following small-p 2 
EM response: 



A x = 0, 



A, 



v 2 



(22) 



Bi 



MS 



[i-(f + f)A 3 ] 
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We see that the image of an insulator is another insulator, 
although with different Hall conductivity and dielectric 
function. 



This transformation, together with eq. (pq), also gener- 
ates an infinite group, denoted Tq(2), which relates dual 
conductors with bosonic charge carriers (such as super- 
conducting films or Josephson junction arrays). 



B. Conductors 



For conductors the form factor III is not analytic as 
p 2 — > 0. The limiting form for small p 2 is related to the 
conductivities by: 



n 3 (p 2 



@xy 



(23) 



p 2 + 



and so the complex quantity II is related to the complex 
conductivity, a — a xy + ia xx , by II = i\fp 2 a* . 

In this case the duality transformation, eq. (jl^), pre- 
serves the momentum dependence of the form factors, im- 
plying II = iyfp 2 a*, with the dual conductivities given 
by the holomorphic relation: 



77 - a 



5 + 77 cr 



(24) 



• Fermion Charge Carriers: If we assume the original 
charge carriers to be fermions, as would be appropriate 
for the integer Quantum Hall systems, then eq. (p4j) re- 
duces to the expression: 



f 



1 -2a' 



(25) 



which is the transformation which has been recognized 
p5[ to imply the experimentally observed duality trans- 
formations p xx — > 1/pxx hi the transitions between the 
a xy — 1 Hall plateaux and the Hall insulator. 

Duality relations between more complicated Hall 
plateaux may be similarly obtained by starting from the 
appropriate anyonic charge carriers, but the same result 
may be obtained more simply by combining eq. ( |25| ) with 
transformation ( fuf ) (which expresses the the absence of 
content of a shift of the statistics parameter by 27r) spe- 
cialized to conductors: 



(26) 



1 - 2a 



Transformations ( |25| ) and (|26|) together generate the 
group, To (2), which is known to relate all of the al- 
lowed (odd-denominator) Hall states to one another, and 
so produces the anyonic Quantum Hall duality relations 
from the simpler fermionic one we have just considered. 

• Boson Charge Carriers: Specializing eq. J24| ) to the 
case 9 = also gives a simple relation between the dual 
conductivities: 



1 

a 



(27) 



C. Duality and Flow 

To this point we have found how the EM response of 
dual systems are related to one another. We wish now 
to understand how the duality transformations change as 
external parameters like magnetic field and temperature 
are continuously varied. The question of how systems 
change as external variables are varied is particularly 
sharp for conductors, since in this case both the origi- 
nal system and its dual are of the same type (i.e. they 
are both conductors). If one imagines a system tracing 
out a curve in the conductivity plane as T and B are var- 
ied, the position of the dual system traces out another 
curve in the same plane. We wish to argue that the re- 
sulting flow commutes with the action of the two duality 
transformations, eqs. (25) and ( |26| ) (or (^6|) and (|2~7j)), in 
the conductivity plane. 

The magnetic field strength, and other other micro- 
scopic properties, enter into the above arguments only by 
changing the values which are taken by the masses and 
other parameters appearing in the effective lagrangian, 
eqs. (|2j) and (^). Temperature partly appears in the 
same way, but also appears in the integration over the 
low-energy degrees of freedom (such as by rotating to 
euclidean signature and imposing periodicity in imagi- 
nary time). Now comes the main point. All of the conse- 
quences of duality follow from the statements = V^ v 
and = P Ml/ , which relate the response functions of 
dual systems (or the analogous statement expressing the 
addition of 2ir statistics flux). And these statements are 
true for all values of magnetic field, temperature etc., so 
long as these quantities do not introduces differences in 
form between the particle and vortex lagrangians, L v and 
r 

Suppose we now consider a system of quasiparticles, 
S(Tq,Bq), chosen for some specific temperature, To, 
and magnetic field, Bq, involving N particles of mass 
m and charge q. Next, suppose that Bq and To are 
chosen to produce the dual system, S, having N vor- 
tices, also with mass m and vortex charge q. That 
is: S(Tq,Bo) = S(To,Bq). Notice that the existence 
of such a To and Bo is plausible given that it involves 
the solution of two equations for two unknowns. For in- 
stance, for particle- vortex interchange the two equations 
are: M(f ,B ) = m(T ,B ) and Q{f ,B ) = q(T ,B ), 
where M and Q are the functions of T and B which define 
the vortex mass and charge. 

Now imagine changing the magnetic field and/or tem- 
perature, say to B' = B (l + S B ) and T = T (l + S T ). 
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The question we ask is this: if we change the dual sys- 
tem by the same amount, to B' = Bq(1 + 5b) and 
T' = Tq(1 + St), is the resulting system still dual to 
the first? That is, is S(f',B') = S(T',B')? The answer 
is 'yes', because the question asked of both systems is the 
same: what is the change in response of a system of N 
objects of mass m and charge q as T and B are varied? 
The only difference between the system and its dual is 
that for S the N objects are particles and for S they are 
vortices. 

Arguing in the same way for the attachment of 2ir 
statistics flux, we see that the entire duality group must 
commute with the flow through the conductivity plane 
as B and T are varied. For fermions this implies a ro(2)- 
invariant flow, while for bosons r#(2)-invariant flow is 
implied. 

D. Beyond Linear Response 

In this section we present a slightly different version 
of part of the previous section's derivation, whose aim 
is to express the action of particle-vortex duality on the 
electromagnetic response without relying on the linear- 
response approximation, eqs. (Q), (|ic|). The argument we 
present assumes fermion charge carriers, since it relies on 
using a statistics parameter 9 = ±7r. Besides clarifying 
our derivation, we present this separate line of argument 
because we believe it will ultimately prove fruitful in its 
own right by explaining the experimental evidence |3r| l 
for current- voltage duality seen in Quantum Hall systems 
beyond the linear Ohm's Law approximation. 

We pick up the story of section (III) just before making 
the linear-response approximations, eqs. (||), ([lO]). The 
response function, T(A), for the system of charged parti- 
cles is given by functionally integrating the fields a M and 
£t weighted by the lagrangian density: 

C(£, a,A) = + ie^OjAoA + C kin (0 + j"(0(<» + 

(28) 

where we have chosen 6 = — n and we need not be con- 
cerned about the detailed form of the particle kinetic 
term, £kin(£)i or current, j^(C)- 

On the other hand, the response function of the system 
related by particle-vortex duality to the original system 
involves integrating the dual lagrangian 

C(y, a, b,A) = - ^ x a^d v a x - e^\d v (a + A) x 

+c kin (y)+r(y)b^, (29) 

where we have instead chosen to represent fermionic 
statistics by choosing 8 = +ir. In order to relate T(A) to 
T(A) we explicitly perform the Gaussian integration over 
a M in eq. ( p9| ) and shift fe M — > + A^ in the result, to get 



C(y, b,A) = + ^ v %d v bx - \ fr x ApduAx 

+C kin (y) + f(y)(b + A) fi . (30) 

To the extent that £kin(y) and £kin(£) (and j^{y) and 
j^(£)) have the same form, eq. (|3(]| ) differs from eq. (|2q) 
only by the term — | & vX A^dyAx, which does not de- 
pend at all on the integration variables. Performing the 
remaining integrations therefore relates T(A) to T(A) by: 

f (A) = T(A) ~ljd 3 x e^A^Ax, (31) 

which is the main result of this section. 

Once specialized to the linear response regime, eq. pi] ) 
implies the relation: 

U!" u {p) = W v (p)+ie^ v p x , (32) 

or: a X y = cr X y + 1. This agrees exactly with what is pre- 
dicted for fermions by a particle-vortex transformation 
followed by a 2ir statistics shift, equation ( |25| ) followed 

V. APPLICATIONS TO CONDUCTORS 

In this section we state some of the observable predic- 
tions which follow from the action of particle- vortex dual- 
ity on the physical systems. We specialize in this section 
to the predictions for two kinds of conducting systems, 
those with fermionic charge carriers (or their images un- 
der repeated duality transformations) - corresponding to 
Quantum Hall systems - and those with bosonic charge 
carriers (or their duality images) - such as for supercon- 
ducting thin films and/or Josephson junction arrays. 

A. Fermions: The Quantum Hall Effect 

We start with Quantum Hall systems, for which the re- 
sults we derive are not new, having been derived from the 
the assumption of duality-invariant flows in ref . pl| . We 
include this case anyway for three reasons. First, this pa- 
per strengthens the theoretical foundation of the assump- 
tion of duality-invariant flow, particularly away from the 
systems' critical points. Second, the experimental success 
of these predictions establishes the existence of systems 
which are clean enough for our quasi-particle/ vortex ef- 
fective theory to apply. Third, we may directly adapt 
many of the Quantum Hall results to the bosonic case, 
for which most of our predictions are new. 

1. Some Group Theoretical Facts 

In previous sections we have found the action of 
particle- vortex duality and 2ir statistics addition to both 
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act on the complex conductivity in a fractional lin- 
ear way, with integer coefficients. That is, we have 
found the transformations to be a subgroup of the group 
PSL(2,Z), defined by: 

a<J + b . . 

cr^ —j, (33) 

with a,b,c and d integers satisfying ad — be = 1. Any 
element of this group can be obtained as products of 
powers of the following two generators: 



S(a) := - i 



T(a) :=a + l. 



(34) 



Figure 1: Flow lines for ro(2)-invariant flow. 

We have found the subgroup of this group which is 
relevant to Quantum Hall systems to be generated by: 



a 



1 - 2(7 

I- a 



= ST 2 S{a), 
= TST 2 S(cr), 



(35) 



1 - 2(7 

or, equivalently, by ST 2 S and T. Any point in the upper 
complex cr-plane can be reached from a fundamental re- 
gion, which we take to be the vertical strip between a = 
and (7 = 1, with the interior of the disc with radius=^ 
centred at a = | removed. As we shall shortly see, the 
boundaries of this fundamental domain are quite gener- 
ally flow lines for ro(2)-invariant flow, and so may be 
seen in Fig. (1). 

It turns out that the subgroup obtained from these 
two generators, denoted To (2), is equivalent to that de- 
fined by requiring the coefficient c in eq. ( |33| ) to be even 
^0| , mjl3| . This alternative definition shows that when 
acting on the real axis, a = er*, it takes rational numbers 
to themselves, with odd-denominator fractions going to 
odd-denominator fractions, and even-denominator frac- 
tions mapping to even-denominator fractions.^ 



"it is tantalizing that the same group, To (2), plays a central 
role in the hierarchical structure of N — 2 supersymmetric 



2. Predictions 

The following consequences follow from the condition 
that a flow in the cr-plane commutes with the group 
To (2): 

1. Universal critical points are predicted for the 
flow at any point, a cr it, which is mapped to itself 
under any member of T (2), "/(a cnt ) = cr cri4 , for 
which the denominator ca + d is neither zero or 
infinite [Tl|Jl^ , ^6| . The complete set of such points 

0~crit — 9 (1 " 



occur at i 



Yang-Mills theories Q. 



r(l+i) and its images under r (2). 
This prediction is borne out experimentally |39|, 
since there is a one-to-one correspondence between 
the experimental critical points of the flow and the 
fixed points of the group. No experimental evidence 
exists for critical points not corresponding to fixed 
points of the group, although these could exist in 
principle. 

2. Although the symmetry does not predict the crit- 
ical exponents at the fixed points [fl9|] , these ex- 
ponents must be the same for all fixed points 
which are related by To (2) pl[^3|Jiq| . This remark- 
able equivalence of critical exponents at different 
fixed points is known as super-universality, and was 
argued for on more model-dependent grounds in 
ref. jh|. Historically, the experimental success of 
this 'prediction' [|d],|^| stimulated the search for an 
underlying symmetry group. 

3. Exact flow lines in the a plane can be derived from 
To (2) invariance plus invariance under particle- hole 
symmetry: a — > 1 — a* . The existence and shape of 
these flow lines depends only on these symmetries 
and not on any other details of the flow's (3 func- 
tion p5| . Figure (1) shows some of the flow lines 
predicted in this way, all of which are semi-circles 
or vertical lines in the cr-plane. The arrows indicate 
the direction of flow to the infrared (whose direc- 
tion does not follow purely on symmetry grounds). 
This very general derivation of the 'semi-circle law' 
- which had been earlier predicted on more model- 
dependent grounds |43| - is spectacularly exhibited 
by experimental systems |ff~4|j . 

4. Flow in the infrared is towards the real axis, ter- 
minating on the real axis at attractive fixed points 
at odd-denominator fractions. Even-denominator 
fractions similarly form repulsive fixed points of the 
flow. One finds in this way a robust explanation of 
the existence of odd denominator Hall plateaux. 

5. Since all allowed transitions between Hall plateaux 
correspond to semi-circles which may be obtained 
by the action of To (2) from the basic semi-circle 
connecting a = to a = 1, one finds a selection 
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rule which expresses which plateaux may be ob- 
tained from which by varying external parameters 
like magnetic field and temperature E3] . The selec- 
tion rule obtained in this way states that a fraction 
P2/92 can be reached by a quantum Hall transition 
from a fraction pi/qi, with qi and qi both odd inte- 
gers, only if \p\q2 —P2qi\ — 1- This agrees precisely 
with all the observed Hall sequences. 

6. There is an element of To (2) which maps each of 
the flow lines to itself, with its ends reversed. For 
instance, the element 7 (cr) = (cr— l)/(2cr— 1), which 
we here identify as the expression of particle- vortex 
duality for fermions, does so for the semi-circle con- 
necting a = to a = 1. When specialized to tran- 
sitions from Laughlin plateaux to the Hall insulator 
this symmetry element is precisely the experimen- 
tally observed |36| duality p xx —> l/p X x with p xy 
fixed H|. 

We regard the natural interpretation of the great ex- 
perimental success of these predictions to be that the 
electromagnetic response of these systems is dominated 
by quasi-particles and vortices. Furthermore the systems 
are sufficiently clean to justify the neglect of those parti- 
cle and vortex interactions which destroy the underlying 
particle-vortex duality, as described here. Finally, the 
relevant quasi-particles are in the fermionic equivalence 
class, in the sense that they are either fermions, or are 
obtainable from fermions by the action of To (2). 

There is also some evidence for a few Hall systems 
for which the critical conductivity is not at the univer- 
sal values flif| , and for direct plateau- insulator transitions 
which do not correspond to semi-circles as predicted here 
fi5f , and an understanding of why particle- vortex duality 
fails here would be very instructive. Since these typically 
involve samples with the most disorder, one possibility is 
that Landau-level mixing is not negligible in these sys- 
tems [ [l6|]47| . We believe this to an instance where in- 
teractions with the disorder ruin particle-vortex duality, 
and so destroy the underlying symmetry of the flow. 

Another potential difficulty often raised in connection 
with this picture is the observed failure of scaling at very 
low temperatures in some samples Q as one passes 
through the critical regime. We put these experiments 
aside, since although these are potentially very telling 
experiments, since scaling is an inevitable consequence 
of a vanishing /3 function, it is not yet clear what their 
proper interpretation is, and indeed there are other ex- 
periments which appear to support scaling, |fl9| J~~| 

There is nevertheless content in the above symmetry 
arguments, since these imply an entire suite of predic- 



'"See, however, ref. ]17| for an alternative interpretation of 
the behaviour near the critical points. 



tions which must all hold together if particle-vortex du- 
ality is valid. So we predict that the above consequences 
of ro(2)-invariant flow should come as a package, with 
the failure of some implying the failure of the others. 

B. Bosons: Superconducting Films 

A fascinating consequence of the generality of the 
particle- vortex duality arguments we present here is that 
they predict different, but equally striking, phenomena 
for the electromagnetic response of other clean two- 
dimensional systems. In this section we describe these 
predictions for systems whose charge-carrying quasi- 
particles have bose statistics (or the image of bose statis- 
tics under a group which we here specify). These pre- 
dictions should have practical applications to supercon- 
ducting thin films and Josephson Junction arrays, and 
some of them have been anticipated flpH for the metal- 
insulating and superconductor-insulator transitions of 
these systems. 

1. More Group Theoretical Facts 

For bosonic charge carriers (and those related to these 
by duality) the action of particle-vortex duality and 
2-7T statistics addition is generated by the following two 
PSL(2, Z) elements: 

a - = ST 2 S{a), 

1 — la 

a -► - - = S{a), (36) 
cr 

or, equivalently, S and T 2 . 

This group is called Tg(2), and is equivalent to the 
condition that a, d are both odd and b, c are both even, or 
vice versa, in the fractional linear transformation, (a a + 
b)/(ca+d). Since bosons may be obtained from fermions 
by shifting their statistics by A9 — 7r, this group may 
be obtained from the Quantum Hall group, Fo(2), by 
conjugating by STS(cr) = cr/(l — cr). Concretely: g £ 
r e (2) implies g = STSh(STS)- 1 , for some h £ T (2). 
(The proof of this statement is easiest to see if the identity 
(ST) 3 — 1 is used.) The simplest way to extract the 
predictions of the group Te(2) is therefore to derive them 
from those of r (2) by conjugating with STS. (Or, since 
S is in Tg(2) anyway, we can equally well get Tg(2) by 
conjugating T (2) with TS(a) = 1- A rather than STS.) 
In particular, it is convenient to choose the fundamental 
region to be the vertical strip between a = and a = — 1 , 
with the interior of the disc with radius=i centred at 
cr = — ~ removed, the boundaries of which again appear 
as particular flow lines in figure (2). 
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Figure 2: Flow lines for rs(2)-invariant flow. Although 
this looks very much like Figure 1, careful examination 
shows different positions for the fixed points, and for the 
directions of flow along the lines. 

2. Predictions 

In this way we obtain the following consequences of 
the commuting of Tg(2) with the flow in the cr-plane, 
which are the direct analogs of those described above for 
Quantum Hall systems. 

1. Universal critical points are predicted for the flow 
at the fixed points of r#(2), which lie at a cr u — i 
and its images, (bd + ac + i) / (c 2 +d 2 ), under rg(2). 
The Ohmic conductivity is always 1 / (odd integer) 
at these critical points while the Hall conductiv- 
ity is even/odd. (These statements are for bosonic 
charge carriers with the same electric charge as an 
electron — these become a cr u — iq 2 and its images 
if the bosonic charge carriers have charge q. In par- 
ticular, the case q = ±2 applies if the bosons are 
Cooper pairs, such as considered in ref. p|.) 

2. The critical exponents at all fixed points related by 
r#(2) must all be the same. In this way the results 
of the scaling theory of ref. j[) for the fixed point at 
o~crit = i may be extended to all of the other fixed 
points which are predicted to exist in the presence 
of magnetic fields. 

3. Exact flow lines in the a plane are immediate conse- 
quences of rg(2) invariance and particle- hole sym- 
metry, independent of the dynamical details of the 
flow's (3 function. These are the images of the flow 
lines of Figure (1) under the conjugation by STS. 
The results are again semi-circles or vertical lines 
in the a plane, implying a new semi-circle law for 
these bosonic systems. (The semi-circle intersect- 
ing the critical point at er = i was anticipated in 
ref. §.) 



Figure (2) shows the flow lines which are predicted 
by the symmetry, with the flow directions given 
which follow from those of Figure (1) (which suc- 
cessfully describe Quantum Hall systems). Notice 
that the resulting flow on the imaginary axis agrees 
with the interpretation of a = i as a metal-insulator 
transition, with a xx increasing or decreasing into 
the infrared on opposite sides of the transition. 

4. More generally, for nonzero magnetic fields the flow 
in the infrared is towards the real axis, terminat- 
ing on the attractive fixed points which are frac- 
tions, a — p/q. The attractive fixed points of the 
flow therefore exhibit the fractional QHE, but with 
fractions for which pq is even (as opposed to hav- 
ing q odd, as was the case for fermions). Fractions 
with odd pq are repulsive fixed points. In partic- 
ular odd integers are repulsive while even integers 
are attractive. 

5. There is a selection rule expressing which fractions 
may be obtained from which by varying external 
parameters like magnetic field and temperature. 
The selection rule obtained in this way states that 
fractions pi j qi can be obtained from p\ / q\ only if: 
(i) pi is odd and q± is even while pi is even and 172 
is odd (or vice versa with the subscripts 1 and 2 
interchanged) and (ii) \p1q2 — Piq\\ = 1- 

6. There is an element of Tg(2) which maps flow lines 
to themselves, with their ends reversed. For flow 
along the imaginary axis the element is simply 
S(a) = —1/(7, or <j xx — > 1/(T XX , which is again the 
expression of particle-vortex duality for bosons. 

We are led to predict the above startling properties for 
bosonic Quantum Hall systems, such as someday might 
be obtained from superconductor-insulator transitions in 
thin films, or from Josephson-junction arrays [^O] 
Unfortunately, these predictions cannot yet be tested 
with ordinary superconductors, because charge carrier 
densities and mobilities in these systems do not put them 
into the quantum Hall regime [p5| . Their verification 
would be the smoking gun for particle- vortex duality, al- 
though its testing must wait until these systems can be 
reliably manufactured in the quantum regime. In this 
way particle-vortex duality might ultimately provide a 
connection between Quantum Hall systems and a xx — > 
l/oxx a $ seen in metal-insulator transitions [5^j5^ ]. 

As in the case of Quantum Hall systems, these pre- 
dictions require very few assumptions beyond the neces- 
sity to be in the quantum regime. First, the electromag- 
netic response must be dominated by quasi-particles and 
vortices, with the quasi-particles being bosons or related 
to these by r#(2) transformations. Second, the systems 
must be sufficiently clean to ensure the absence of interac- 
tions which distinguish the quasi-particles from vortices, 
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and so ruin their similarity. Subject to these conditions 
all of the above consequences follow as a package from 
the quasi-particle- vortex duality symmetry of the flow in 
the conductivity plane. 

VI. CONCLUSIONS AND OUTLOOK 

We have argued that a broad class of dual relationships 
arise in two dimensional systems for which the EM re- 
sponse is governed by particles and vortices whose prop- 
erties are similar (perhaps because they are weakly in- 
teracting). For systems having fermions as the particles 
(or those related to fermions by the duality) the particle- 
vortex duality implies the duality group is a level-two 
subgroup of PSL(2, Z) called r (2). We argue that this 
duality has been observed, since this group has been 
previously identified as explaining many of the observed 
properties of the EM response in Quantum Hall systems. 

The generality of our arguments lead us to propose 
a similar, but distinct, set of properties for clean two- 
dimensional systems whose EM response is controlled by 
bosons as the particles (or those related to bosons by the 
duality) . The group implied in this case is another level- 
two subgroup of PSX(2, Z) called Fg(2). The observation 
of these specific predictions in these systems would be the 
definitive test of our ideas. 

One might ask whether other duality symmetries apart 
from those described above might arise in other systems. 
More concretely, there are precisely five level-two sub- 
groups of SL(2,Z) ]lOfl , so one might wonder if other 
choices for quasi-particle statistics might generate these 
other three groups not yet used J5f|. In fact, two of these 
groups can be ruled out as symmetries acting on the com- 
plex conductivity, because they do not contain the gen- 
erator corresponding to the addition of 2ir statistics flux. 
However, because these groups are related to To (2) and 
rg(2) by conjugating by S, they can be thought of as the 
action of these latter groups in the complex resistivity, 
rather than conductivity, plane. The third group, T(2), 
is contained in the other four, and has been proposed 
elsewhere to play a role for Quantum Hall systems [j59| , 
in particular when the splitting between electron spins 
is much smaller than the gap between successive Landau 
levels @. 
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APPENDIX A 
The Statistical Gauge Field 

In this appendix we review the concept of the statisti- 
cal gauge field p5[ . This serves not only to remind the 
reader of the construction, but also to set up the nota- 
tion. This appendix is based on the review article by 
S-C. Zhang, || . 

• First- Quantized Formulation: 

Consider a configuration of N charged particles, each 
with the same charge e and mass m moving in a two di- 
mensional plane, with positions xi, . . . , x^r and separa- 
tions Tij = Xi — Xj = —Yji. Allowing for particle-particle 
interactions, e 2 V(~x.i — x^), and the possibility of random 
static impurities giving, rise to an electric potential U (r), 
the Hamiltonian for the system can be written as: 

H = 2^ £ ir im - - eA *^) 2 + e2 *,•) 

i 

+e^C/(x t ), 

i 

where a = 1,2 and we use units in which the speed of 
light c = 1. One can include a neutralising background 
field if desired, without changing any of the subsequent 
analysis significantly. 

Following MacDonald j6l]], a gauge transformation to 
a new Hamiltonian is defined as follows. Let tp^j be the 
angle between the vector and an arbitrary fixed di- 
rection, e.g. the x-axis, so that ifij — > (pij + n if the two 
particles i and j are interchanged. A gauge transforma- 
tion from the multi-particle Schrodingcr wave-function, 
tpix-i, . . . , x„) to a new wave- function ip(xi, ■ ■ ■ > x «) is 
defined as follows, 

V>(xi, . . . , X jv) =e TZ ^<i' Py V(xi,...,Xiv) (38) 

for a constant 9, as yet arbitrary. Under interchange of 
any two particles, i and j, the phase factor changes by 
e l8 . Thus, if 9 = 2kn for some integer k, the phase factor 
is unity and the new wave function has exactly the same 
statistics as the old one — if it was symmetric under 
interchange of two particles before it remains so, if it was 
anti-symmetric under interchange of two particles before 
it remains so. If 9 were an odd multiple of n a bosonic 
wave function would be transformed into a fermionic one 
and vice versa. 

The gauge transformation ([38]) can be incorporated 
into the vector potential A^, fx = 0, 1, 2, for the electro- 
magnetic field by defining a new field (the statistical 
gauge field) 
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h 9 



a Q ( Xl ):= V V^' 

e 7r * — ' 



ifij 



(39) 



e af, (doa l3 -dpao) = -=^j a 



2W 



(46) 



where the gradient operator acts on the position of the 
i-th particle. The new gauge transformed Hamiltonian is 
then 



This equation can be obtained, together with equation 
@), by treating dynamical field in the action 



H 



-E (- 

2m 4^ V 



5* = 



eA a (-*Li) - ea Q (x j; ) 



5 ^ , hH , , dtd 2 xe^ x a^d v a x + \ dtd 2 xa^, (47) 



+e 2 Ei<i V(x t - x,) + e Ei U{x<). (40) 

It is stressed that the physics of this new Hamiltonian is 
identical to the old one, provided is an integer. 

• Second- Quantized Formulation: 

We now reformulate the above first quantised treat- 
ment in second quantised form, ignoring spin. Accord- 
ingly define, 

H = 

/ d 2 x*t( x ) \^(-ih V - eA(x) - ea(x)) 2 + e£/(x)l (x) 
+ | / d 2 xd 2 x'5p(x)V(x - x')fy(x'), 

where p(x) = e^{x)^f(x) is the charge density and 
5p = p— < p >. 
Equation ( |39| ) can be expressed as 



which is the Chern-Simons action for the statistical gauge 
field coupled to a source. 

The problem can now be re-written in path integral 
form. Define the original Lagrangian density for the mat- 
ter fields 



= * f (ihd t - eA 



— i-ihV -eA) 2 -eU 
2m y — ' 



1 



d 2 x'5p(x.)V(x - x')fy(x'). (48) 



After introducing the statistical gauge field the new La- 
(41grangian is 



= &{iW t - e(A + a ))# 
1 



2m 



eU 



a„(xf) = 
which in continuum form is 



?i0\ ^ (r 

/ tap — 



(42) 



ih V — eA — ea)^ 
~[d 2 x'6p(x)V(x-x')6p(x% (49) 
together with the Chern-Simons Lagrangian 



a Q x) = - — 
;From this follows 



d 2 x> e -4^ 



-p(x'). (43) 



Co 



(50) 



Then the effective action for the gauge field is ob- 
tained, in the original formalism, from the path integral 



e^V^x) = ^p(x), 



(44) 



(51) 



since V 2 ln(|x - x'|) = 2tt^ 2 )(|x - x'|). The curl of a 
is non-zero here, despite it's original definition as a gra- 
dient, because of Aharanov-Bohm type singularities in 
(|39|). These manifest themselves as (5-function singulari- 
ties in the curl of a in the first quantised theory but are 
spread into a smooth distribution on the right hand side 
of , in the second quantised theory. Equation (|44|) is 
a dynamical equation of motion for the statistical gauge 
field a a, but it is not yet co- variant, as it lacks a terms. 
To include these allow a a and p to depend on time, then 



differentiating (44) and expressing the result co-variantly 
gives 



e p Ofj(d a a - d a a Q ) = —p = T d a j a , 

where j a is a current. This integrates to 



(45) 



where 5^ = f dtd 2 x£y. On the other hand, after intro- 
ducing the statistical gauge field, the partition function 
is 

Z[A] = J Va li V^V^e iSa[a ^ +iS ^ A " +a ''^ t ' w] 

= J Va^ iSa{a ^Z[A + a], (52) 

where S a = J dtd 2 xL a . 

So far all manipulations have been exact. For 6 — 2kir, 
with integral k, the theory described by (|5^) is identical 
to that with no statistical gauge field at all, i.e. Z[A] — 
AfZ[A], where N is an irrelevant constant that will be 
ignored in the sequel. 
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APPENDIX B 
The Relativistic Trick 

We now follow Kivelson, Lee and Zhang jl2) and in- 
tegrate out the matter fields. This cannot be done ex- 
actly, of course, but on general grounds one expects the 
effective action to be invariant under local gauge trans- 
formations. Firstly ignore the statistical gauge field and 
consider the partition function (|o"l|). Naive power count- 
ing arguments then imply that the most relevant terms, 
at least in the infra-red limit, are 



FefM] 



d 3 xd 3 x' 



Fj (2)71-1(2; - x')Fio(x') 



F 12 (x)n 2 {x - x')F 12 (x') 



\e^ x A^x)^{x-x')d v A x {x') 



(53) 



where Fio = Ei, i = 1, 2, is the electric field and F12 = B 
the magnetic field (the measure here d 3 x is a shorthand 
for dtd 2 x and Greek indices, fj,, v, A take three values, 0, 1 
and 2). Note the presence of the non-local form factors 
7Ti , 7T2 and 7T3 , where the argument x — x' includes time 
as well as space. 

There is an implicit assumption here, that the terms 
that are most relevant by naive power counting are also 
the most relevant in the full theory, at least in the infra- 
red limit. This is a very strong assumption, as naive 
power counting gives the most relevant operators of the 
free field theory and we have here a strongly interacting 
theory. One of the underlying assumptions of Kivelson, 
Lee and Zhang's argument is therefore that any anoma- 
lous dimensions in the strongly interacting theory do not 
change the conclusions of naive power counting. In fact 
the form ( |53"| ) is rather more general than that — it will 
be argued below that this is the most general possible 
form in a momentum expansion, independently of any 
power counting arguments, provided one allows the po- 
larizations to depend on the magnitudes of the external 
fields, E 2 and B 2 . 

The mathematical manipulations involving equation 
( p3| ) can be simplified by writing it in relativistic form 
using the following trick. We introduce a "metric" on 
3-dimensional space-time defined by 



9nv ( x 




(54) 



and write 



d 3 xd 3 x' I - x r )^-det{g)g^g^F^{x)F pT {x') 



-\e pvX A p {x)^(x~x')d u A x {x') 



(55) 



where all metric components are functions of x — x' . Note 
that the metric does not appear in the Chern-Simons 
term at all — as is well known it is independent of the 
metric. 

The calculation can be simplified by setting 7Ti = 7:2 
and working in the relativistic formalism with the 3- 
dimensional Minkowski metric, r] pl , = diag(—l, 1, 1), and 
relativistic action 

reff[A] = 



d 3 xd 3 x' { - - 7ri(z - x')F^{x)F liV {x l ) 



^ vX A^x)n 3 (x - x')d v A x (x')\ , (56) 



where F pv = rf p r\ vr F pT . Provided everything is kept co- 
variant the non-relativistic expressions can be recovered 
by re-instating the metric ( |5~4| ) at the end. From now on 
we shall use the simpler relativistic notation of equation 

In relativistic notation we can argue that ( |56"| ) actu- 
ally encodes higher order terms too. In strong fields 
one would expect terms such as F p V F V p F p ' X F X p , and 
other Lorentz invariant powers, to be present in the effec- 
tive action even in the low momentum regime — though 
terms like d 2 F^ v F llv arc definitely ignored. In three di- 



mensions terms like F p V F V p F p X F 



x p can be absorbed 



into ( |56| ) by allowing the polarizations to depend on the 
Lorentz scalar F pv F pv . This is because we can always 
exchange the anti-symmetric tensor F pu for the vector 
F p = ( pvp F vp and the only way to make a Lorentz scalar 
from products of F Pl ■ ■ ■ F^" is to contract the indices in 
pairs (and n must be even) , so all such terms can simply 
be incorporated into the statement that tt\ depends an- 
alytically on the Lorentz scalar F 2 = F pv F pv and then 
Taylor expanding iri in powers of F 2 . An exactly 
similar statement applies to ^3. We shall continue with 
the form ([56]), bearing in mind that when the external 
fields are strong the polarizations may depend on them 
explicitly. 

It will be more convenient to work in momentum space 
where (pq) reads 



r eff [A) = jd 3 p{-[- A ) w^F^i-p^ip) 

-^ x A^-p)iT 3 {p)p u A x {p) \ , (57) 



with d 3 p 



d 3 P 



"In general this would necessitate the introduction of multi- 
point interactions with n points, Xi, . . . ,x„, but in the long 
wavelength limit two point interactions suffice to extract 
conductivities. 
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A more compact way of writing fl57|) is 



r eff [A]=-J d^jillj 



(p)A^-p)A^A u (p) 



+-IL s (p)J^A^-p)A u {p) } , (58) 



where A% = 5% - J» v = ie^px/yjp 2 , ^ = 

Ttip 2 /h 2 and II3 = nzy^pP/h. This form is useful, 
since A'< Q A% = A^, A^ a J au = J^ a A a u = J» v and 
J^ a J av = — A M „, which makes Gaussian integrals par- 
ticularly simple. Note that J is a Hermitian matrix for 
space-like momenta, but anti-Hcrmitian for timedike mo- 
menta — however \[p 1 ~ J is always Hermitian, which is all 
that is necessary to ensure that the momentum space ex- 
pression is always Hermitian. 

It is stressed that (^) is the most general form of 
the effective action in the long- wavelength limit, even for 
strong fields, provided Hi and n3 are allowed to depend 
on |E| 2 and |B| 2 . 

APPENDIX C 
The Duality Prescription 

The purpose of this appendix is to derive the dual- 
ity transformation from eq. (||) to eq. (|4|), which we do 
following the general duality prescription |?0| . Our start- 
ing point is an extended Lagrangian which is obtained by 
coupling to a new gauge field, A^, which is constrained 
to be pure gauge: 



— e"** ail d x a v + £ p (t,a + A) 



20 

K 



[d^-q^ + A^+A^Y 
+e^%d„A x + ■■■, 



(59) 



where the field b^ is a Lagrange-multiplier field which is 
introduced to enforce the vanishing of the field strength, 
d^Av di/Afj, . 

That this extended action is precisely equivalent to 
the original action may be seen by first performing the 
integration over 6 M , which produces a functional delta 
function which enforces the constraint e^ vX d u A\ = 0. 
This, together with the gauge fixing condition, implies 
that the integration over A^ is equivalent to setting _4 M — 
everywhere in the path integral, which reduces eq. ( |59| ) 
to (||), as claimed. 

The dual version of the lagrangian is found by instead 
performing the functional integrals in a different order, 
integrating out cj) and A^ and leaving as the dual vari- 
able. Care must be taken when performing these inte- 
grals to properly handle the vortex boundary condition 
which is satisfied by 4>, namely: 



0(0 + 2tt) = 0(0) + 2vr q4> ^ N a , 



(60) 



where N a are integers labelling the vorticity of each vor- 
tex, and $ is the angular polar coordinate taken at spatial 
infinity, a long distance away from the vortex positions. 

To integrate over (j) it is convenient to write (f> = lj + (p, 
where w is a particular configuration having the same 
boundary condition as does </>, so ip is simply periodic: 
ip^ + 2tt) — <^(f?). For u> we choose 



w(x) 



2tt 



E 



N a arctan 



2 

Va 



(61) 



where y l a ,i = 1,2 are the coordinates of the positions 
of each vortex, with the index i labelling the two space 
directions. 

Notice that the gauge potential defined by t> M = d^uj, 
has vanishing field strength, except at the positions of the 
vortices, where it has <5-function singularities, so: 

e^b^vx = -—J2 N a y^{t)b^[x - y a (t)}. (62) 



This term appears in the dual lagrangian, and provides 
the minimal coupling of the vortex positions to the po- 
tential & M . 

With these definitions, the integrations over ip and A^ 
are straightforward. It is simplest to choose ip = as 
a gauge condition, and then directly perform the uncon- 
strained Gaussian integral over A^. The result is eq. (Q), 
without the vortex kinetic term. 



APPENDIX D 
Particle- Vortex Duality and Landau Level 
Addition 

In this appendix we derive the duality transforma- 
tion between pseudo-particles and vortices, but within 
a second-quantized path-integral framework. (See also 
ref. |(|| for a more detailed discussion of particle/vortex 
physics within the abelian Higgs model supplemented by 
a Chern Simons term.) The treatment here is based on 
f32f , except that we use a relativistic notation since this 
makes the manipulations simpler and, as shown in Ap- 
pendix B, the non-relativistic form is easily recovered 
from the relativistic form. 

We start from the action for a complex scalar field $, 
with charge e, coupled to an external electromagnetic 
field Afj_ with a statistical gauge field 



S[$,A,a] = 



with 



d 3 x 



29h H v p 



S m [$,d] (63) 



-\Jd 3 x [(i^ - ea M )$] f [(ihd" - eo")$] + S mt [\<S>\ 2 } (64) 
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where a = A + a, and Si n t is an interaction term, possibly 
including a mass term. 

• The Flux Attachment Transformation: 

The effective action involving the statistical gauge field 
should reproduce the same physics as (|58|), if 6 = 2kir. To 
examine this further we first observe that the statistical 
gauge field a only ever appears in the action (|64|) in the 
combination A + a, so define a = rjA + a , (at the moment 
r/ = 1, but it is introduced here for later convenience). 
Now integrating the matter fields out of ( |52| ) gives rise 
to an effective action for A^ and a M of the form 



r$ f [A,a] = - / d A p\ |n 1 ( P )o / ,(-p)A'"'a I/ (p) 



i. 



(in. 

\2h9 



--n 3 ( P )J' u 'a M (-p)o v (p) 



d 3 p^-J^a^-p)a v (p). (65) 



Even though the matter integrations cannot be done ex- 
plicitly, the construction ensures that the form factors 
III and II3 appearing in equation ( |65[ ) must be identical 
to those appearing in fl58|). They will be functions of the 
field strengths for a^ in general but, provided a M is small 
(i.e. provided a M almost exactly cancels the external field 



rjA^), they can be evaluated at 



0. 



Now equation (|65|) should describe exactly the same 
physics as (58), after a is integrated out, and this is what 
gives rise to the flux attachment transformation, as we 
now describe. 



Integrating the statistical gauge field out of (65) is eas- 
ily achieved when a M is set to zero in III and II3, as it is 
then quadratic in a M . The resulting effective action for 
A^ will have the same form as ( |58|) with III and II3 eval- 
uated at A^ — 0, but with different form factors, which 
we shall denote by II 1 and II3. 

A subtlety in the integration is that in the quadratic 
from a^M^ a v the matrix M pv is not invertible as it has 
a zero eigenvalue. However it is only really necessary to 
find a matrix M. such that, in relativistic formalism in 
momentum space, M. jlv M vp = 5^ p — ^p- in order to 
integrate out a. The result is 



Veff[A] 



d*v\~ni(j>)A^-p)k^A v (p) 
+ lfl 3 (p)j^A^p)A u ( P ) 



(66) 



where 

IT = 

fia = 



e 2 7]Tr\ 2 p 2 Ui 



h9 



h 2 D 



/ /II 

1 jr / e T) 7r 



he 



e 2 r]ir 

he 



n 3 + ^(^f) 



D 



(67) 



with 



D 



n 2 



El 

h 



n. 



p 2 2 

l-Trf + I 7T.3 



h6 



he 



(68) 



The non-relativistic transformation is obtained by 
re-introducing the metric ( |54| ) and noting that p 2 = 

become 



-Po 



f^p , and then equations (|£ 



7T2 



71"3 



e 2 r]ir^ 


2 

7Tl 


he ) 


d 


e 2 r/TT^ 


2 

7T2 


he ) 


d 


e 2 r] 2 Tr 


) 


he 





e 2 ?77r 
/i0 



(69) 



where d = (—irfp 2 , + ni7r 2 p 2 )/h + ( 7r 3 + %g L ) , and this 
is the form of the transformation given in (I4] (except it 
is given in Euclidean signature in that reference). 

For a Hall conductor the transverse conductivity is re- 
lated to the polarization tensor n3 by 



Cxy — 7T3 = MI3/1 



(70) 



while the transverse conductivity involves breaking the 
3-momentum up into frequency and spatial momentum 
= (hu, p) and taking the limit 



= lim 



' p 2 TT\/% 



= lim 

p=0 



p=0 



= -^im [ni/cj] p=0 . 



(71) 



The transformed conductivities, a xx and cr^ are related 
to the transformed polarization tensors Hi and n3 in a 
similar way. 

Equation ( |69| ) then gives, with 77 = 1, 6> = 2fc7r and 
units in which e 2 //i = 1, 



4fc 2 (o- 



2 



• o" 2 3/) - 4ka Xi 



2k(a 2 xx + a 2 ) + o 



xy 



^k 2 {cj 2 xx + a 2 )-Akcr xy + l 



(72) 



It is stressed here that the above analysis represents a 
symmetry under certain conditions, as discussed in | f[2j , 
such as very low temperatures. In the long- wavelength, 
zero frequency limit, the phase diagram of the quantum 
Hall effect is symmetric at low temperatures under the 
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above transformation. For example the critical exponents 
at related second order phase transitions should be iden- 
tical. 

It should be remembered that a xx and a xy represent 
the components of a tensor, it may at first sight seem 
unnatural to be applying a non-linear map which mixes 
up the different components of a tensor — what about 
co-variance of the tensor components? In fact equations 
( f72| ) are very natural from this point of view. If we define 
a complex co-ordinate z :— x + iy, and it's conjugate 
z = x — iy, then the conductivity tensor in these co- 
ordinates is reduced to a single quantity, a :— a xy + 
iffxx, with a positive imaginary part, since a xx > 0. The 
transformation reduces, in this co-ordinate system, to 



a = r) 



1 + ^ 



(73) 



v(r) 



(77) 



where (ya,y 2 ) denote the position of the vortex labelled 
by a, which should be summed over in the path integral. 
In these variables the matter action reads 



S m [$ ,&,a] = -- I <l\v 



$o(^<9 M i9 + ifiv*d ll v - ea M ) 



s int \m 



where indices are understood to be contracted with the 
Minkowski metric (— 1, +1, +1). 

We can now perform the integral over i? by introducing 
an auxiliary vector field, A M — Tid^d, and imposing the 
constraint d^Xu — d^X^ — with a Lagrange multiplier 
field, bp. So we write 



which gives a = 1 ^ 2ka for r\ — 1 and 9 = 2kn. This last 
form can be obtained by k iterations of the generating 
transformation 



(74) 



l + 2er 



as is easily checked. Equation (74) is the transformation 
ST- 2 S in the text. 

• Particle- Vortex Duality: 

We could consider eq. (^34j) as either: i) a bosonic prob- 
lem which is transformed from another bosonic problem 
(9 = 2kir) or ii) a bosonic problem transformed from a 
fcrmionic problem (9 — (2k + l)7r) — only in the for- 
mer case is the statistical gauge field transformation a 
symmetry. 

It will be argued in this section, following [pjj , that 
there is a second symmetry in the bosonic case (fc even), 
a — > —1/<t. This is a Z2 symmetry which maps an 
insulator a = to a superconductor a — ioo and has 
o~ xx = e 2 /h, <j xy — as a fixed point. 

We derive the second duality transformation as follows: 
firstly write (63) in terms of the pseudo-particle para- 
magnetic current, j M = - (9,,$)+$], as 



S[$,A,a\ = - J d 3 



20h 



dfx 



-,2 



+ S int [\$\ 2 ]. 



(75) 



Alternatively the action ( p3[ ) can be written in terms 
of the vortex current by splitting $ up into a smooth part 
and a vortex part as 



$(r) = $ (r)e- l,J(r) w(r) 



(76) 



where $o( r ) is real, $(r) is real, positive and single valued 
and 



Performing the functional integral over A M puts ( |63| ) in 
the form 



S[Q ,A,a,b] = / d 



- — — & vp a d a 
29 h °" v p 



1 



~CCf b r 

sUHl (so) 



where = dpb v — d v b^ is the field strength for 6 M and 
the vortex current is = ^e^ x (d v v*)(dxv). The in- 
tegration over A also induces a functional determinant, 
ln(det$o) and this has been absorbed into the interaction 
term for $0 , as indicated by the prime on S' int [$q] . 

Now integrate a M out of (p0|). The only terms involving 
a are 



Sa 



d 6 x 



29 h 

So integrating out a we get a term 
h 2 



d 6 x 



2n e 



(81) 



(82) 



The action in terms of b is now 



SM[$ ,A,b] = Jd 3 



2ne 2 » up 



d 3 x 



4e 2 $2^ 



- / d 3 2 



(83) 
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Now let b = b + jr^A and the action becomes 



d 3 x 
d 3 x 
d 3 x 
d 3 x 



7r e 



e^b^bp 



29 h 

7T h 

29 e 2 

2 



5Jnt[*o]. (84) 



where 6* = — 

For comparison ( f75| ) reads, writing $ = ^oe - '" 5 so 
that the pseudo-particle para-magnetic current is j M = 



5 (a) [$o,A' 



cfai 



29 h ^ p . 



e 2 $n 



The duality symmetry that we seek lies in the sym- 
metry between (H) and (|85|). One way to see this is to 
fix the external field A^ at some background value and 
treat <&o as a classical field. If $o is non-zero constant 
the action ( |S5|) describes a superconductor while (Q) de- 
scribes an insulator — as argued in fjj vortices condense 
to give an insulator. If the gauge symmetry is not broken 
and describe conductors dual to one another. 
To examine this situation we analytically continue $q to 

complex values and set e 2 <i>g w ^-^-C = -r-wC where 
C is dimensionless (this is the leading term in the Fourier 
expansion of <I>q, in the long wavelength limit when the 
gauge symmetry is not broken). So the longitudinal con- 
ductivity is (setting p^ = (twj, 0)) 



n2 — 



lim 



(86) 



(see equation @ with tt 1 = fi 2 e 2 $l/p 2 ). Similarly S 1 ® 
describes a conductor with longitudinal conductivity 



lim 



1 h 

C^ 2 



(87) 



(equation Q) with m = l/(e<J> ) 2 ). 

For a fixed external field <&o, and so the conductivities 
a xx and <J XX , depend on the external field as well as other 
external parameters, such as the temperature. In general, 
therefore, it should be possible to find pairs of values for 
the external parameters, labelled generically by X, such 
that C(X) — l/C(X'). (this is equivalent to setting 
V = P in section IV). Then the effective actions (p3) 



and ( j35| ) are identical except for the extra Chern-Simons 
term for A^ in 

In particular if we take 9 = —n in J8^| ) and n in ( |85| ) 
(these two values of 9 are of course indistinguishable in 
the statistical gauge field transformation) the sole effect 
of the extra term on the effective action is to shift the 
transverse conductivity by a xy — > o xy + 1 (where we have 
set e 2 /h — 1). Of course 9 = ±ir means that we started 
from fermions and transformed to bosons using the statis- 
tical gauge field so we can argue that we have derived the 
T transformation (which corresponds to Landau-Level 
addition in the first quantised theory) as follows: start 
with a fermionic system and transform to a bosonic sys- 
tem using 9 = ±7r, then use the above argument to show 
that tr — ► (7+1 is a symmetry and then transform back to 
fermions. This shows that o — > a + 1 is a symmetry for 
the fermionic system (it is not a symmetry for a bosonic 
system, because 9 — ±ir does not keep bosons as bosons). 
This argument, while plausible, has ignored the currents 
and the functional integration over <!>o. 

A more convincing argument takes the currents, both 
pseudo-particle and vortex, into account. Pseudo- 
(8$>^irticles and vortices are massive and so should decouple 
in the long wavelength limit. Just as in the derivation of 
the flux attachment transformation we can argue that, in 
the long wavelength limit, integrating out $0 from fl85| ) 
and summing over pseudo-particle currents must lead to 
an effective action for the external gauge field of the form 



d 3 xd 3 x' ( (\ ) (/<V(zK(x - z'X/TV) 



+ ~e» v %(x)Tr%(x-x')d v a x (x') 



witri fjlv = d^d v - d u a^ = d^(a v + A v ) - d v (a^, + A M ) 
(we have set tt% = 7r" in the relativistic form — the su- 
perscript a indicates that these are polarization tensors 
associated with the field a M ). As before integrating out 
a now gives an effective action for the external field A, 



Teff[A] 



d 3 p||aif 1 ( P )^(-p)A'"'^(p) 

+ ^n 3 (p)J^A IM (-p)A v (p)\, (89) 



with polarizations 



713 (w,p) 



/ e 2 7r^ 


r- 1 




' d a 


e 7r 


/e 2 7r 


~h9 ~~ 


\ h9 



+ (—) 
X \ he > 

d a 



(90) 
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where d a 



{■Ktn-Pl + P 2 )/h 2 + (tT§ + (^)) 2 . This 



leads to complex conductivities that are related by 

a a 



(i + fO- 



with a a = lim w ^o(7T3 — oJirf) and e 2 /h = 1 (remember 

lim w ^ (w7T?) = -i<?Z x )- 

Similarly integrating out $o and summing over vortex 
configurations in ( |84| ) must lead, in the long wavelength 
limit, to 



r ( 3 f [A,b} = 



d 3 xd 3 x' 



d'x^PA^Ap 



7T h 

26 e 2 



d 3 xe^%(x)dM(x). 



(92) 



Remember that here = 9^(6,, + ^-fA v ) - d u (b^ L + 

sin. a ) 

Now we can integrate b out of ( |92| ) in exactly the 
same way as a was integrated out of ( |88| ) — the only 
differences are i) the presence of the term involving 
— ^^& vp A^dyAp which subtracts from the right 
hand side of 7r 3 in (|9p|), ii) ^ is replaced with ^ with 

= — 2? in (90); iii) we should use 77 = ^ and iv) the 
7r" are replaced by ir b , which will in general be different. 
The resulting effective action is 



t e ff\A) = -\ J d^p^^A^-p^Auip) 

+^-h(p)J^A^-p)A„( P )\, (93) 



with different polarization transformations to (| 
equation (p9h, 



see 



, 9 / \ 2 b 



7ri(w,p) 



2 \ 2 



d b 



9 

e 7r 
(94) 



where d b = (7rJ) 2 (-pg + p 2 )/?i 2 + (tt| - (^)) 2 . 

In terms of complex conductivities ( j94[ ) results in 
(again with e 2 //i = 1) 



(95) 



i° h --Y 



with cr fc = lim w ^o( 7r 



If 6> 



27rfc then this 



equation represents a symmetry of a bosonic theory, 
a — > — —z^k w hich is equivalent to the modular trans- 
formation 7 = ( ^ 2k) P ar ^^ cmar case fc = is 

As discussed in Appendix B the effective actions ( |89| ) 
and (|95| ) are long-wavelength but not necessarily weak 
field if the polarisations depend on the background elec- 
tric and magnetic fields. This means that the conductivi- 
ties a a and <r b will depend on external parameters, gener- 
ically labelled X before, including the external magnetic 
field and the temperature, so consider ([)3|) to be the ef- 
fective action for a perturbation on the external field. 
Typically it will be possible to find pairs of parameters 
so that o a (X) = a b (X'). Then eliminating a a and a b 
from @) and @ gives S (a = - 1 /a) when 6 = 0. To- 
gether with the flux attachment transformation (|74|) this 
generates the group Tg for bosons. Conjugating with ST 
as described in the text gives To (2) for fermions. 

APPENDIX E: 
The Gaussian Integrals 

In this appendix we perform in detail the Gaussian 
integrations which lead to our central formula, eq. (fl5|). 
Our starting point is eq. (|J) after the integration over 
the positions of the quasi-particles and vortices has been 
performed, using eqs. (|) and (Q. This starting point 
may be written: 



d 6 x 



26 



^a^xau + e^b^a + A^ 



- I d 3 xd 3 y (a + A) fi (x)P^(x-y)(a + A) u (y) 



+b fl (x)V^(x-y)b u (y) 



(96) 



Our task is to perform the integrals over a p and 6 M , 
and our interest is in the dependence of the result on A^. 
Because none of the functional determinants which arise 
in the integrations depend on A^, these may be neglected 
as contributing an ^^-independent additive constant to 
the EM response function, T[A]. 

Performing first the & M integration, the saddle-point 
condition may be written as: 

b^x) = - [d 3 y V^(x - y y af3 d a (a + A)p{y), (97) 



where the kernel satisfies the definition 

d 3 y V^(x - y)V x »{y - z) = A/ (x - z), (98) 



for the projection operator 

A/ Or -z) = (S; ~ d^/d 2 ) 5 3 (x - z). 

The result for the effective action then becomes: 
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Si=- 



where: 



(99) 



d 3 xd 3 y (a + A) u (x)P^(x -y)(a + A) v (y), 



ot@v 



dp. 



(100) 



The construction of kernels like V^ v is dramatically 
simplified by working in momentum space, and expand- 
ing the kernels in terms of the basis of tensors 
and J^ v , as introduced in the text. For instance, if 
V^(p) = AxW+A 2 J» v , then V uv (p) = B 1 K^ + B 2 J^ V 
with 



Bi 



A\ + AV 



Bo 



A 2 



Aj + Ai 



(101) 



which is more compactly written in terms of the complex 
variables A = Ai + iA 2 and B = B\ + iB 2 as: B = 1/A. 
Proceeding in the same vein, the saddle point for the 
integration is 



d 3 yK,^P vX A x , 



(102) 



where K^K XlJ = A/ and K» v = P*" + % e» Xv d x . Using 
this in the integration over a M gives the electromagnetic 
response function of eq. ([ll]) , with the EM response ten- 
sor, given by: 



IF" = P*" - P» x K,xpP p " 



(103) 



Once expressed in terms of complex variables, as 
above, this result gives eq. (|l5|) of the text. 
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